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Abstract. 

X-ray satellites since Einstein have empirically established that the X-ray lumi- 
nosity from single 0-stars scales linearly with bolometric luminosity, Lx ~ 10~^ L^oi. 
But straightforward forms of the most favored model, in which X-rays arise from 
instability-generated shocks embedded in the stellar wind, predict a steeper scaling, 
either with mass loss rate - M ^ L^J^ if the shocks are radiative, or with ~ 
~ L^^^i if they are adiabatic. We present here a generalized formalism that bridges 
these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to 
the local radius. Noting that the thin- shell instability of radiative shocks should lead to 
extensive mixing of hot and cool material, we then propose that the associated soften- 
ing and weakening of the X-ray emission can be parameterized by the cooling length 
ratio raised to a power m, the "mixing exponent." For physically reasonable values 
m ^ 0.4, this leads to an X-ray luminosity Lx ~ M^'^ ~ L^oi that matches the empirical 
scaling. We conclude by noting that such thin- shell mixing may also be important for 
X-rays from colliding wind binaries, and that future numerical simulation studies will 
be needed to test this thin-shell mixing ansatz for X-ray emission. 



1. Introduction 

Since the 1970's X-ray satellite missions like Einstein, Rosat, and most recently Chan- 
dra and XMM-Newton have found hot, luminous, O-type stars to be sources of soft 
(< 1 keV) X-rays, with a roughly linear scaling between the X-ray luminosity and 
the stellar bolometric luminosity, L^c ~ 10~^ L^oi (iGtidel & Nazell2009l ). In some sys- 



tems with harder (a few keV) spectra and/or higher Ljc, the observed X-rays have been 
associated with shock emission in colliding wind binary (CWB) systems, or with mag- 
netically confined wind shocks (MCWS) (see reviews by Corcoran and Gagne). But 
in putatively single, non-magnetic O-stars, the most favored model is that the X-rays 
are emitted from Embedded Wind Shocks (EWS) that form from the strong, intrinsic 
instability (the "Line-Deshadowing Instability" or LDI) associated with the driving of 
these winds by line- scattering of the star's radiative flux (see review by Sundqvist). 

This LDI can be simply viewed as causing some small (< 10"^) fraction of the 
wind material to pass through an X-ray emitting EWS, suggesting then that the X-ray 
luminosity should scale with the wind mass loss rate, ~ M. But within the standard 



Castor et alJ (119751 . hereafter CAK) model for such winds, this mass loss rate increases 
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Figure 1. Illustration of cooling zone within a wind shock, showing associated 
scalings for X-ray luminosity with mass loss rate M and bolometric luminosity 
Lbou for adiabatic shocks with cooling length much larger than the local radius, i » 
r, or radiative shocks with £ <^ r. Thin- shell mixing of such radiative shocks is 
posited to lead to a reduction of the X-ray emitting fraction that scales as a power- 
law of cooling length, fx ~ f'^. For CAK wind index a, a mixing exponent m = 
I - a leads to the observationally inferred linear scaling of X-rays with bolometric 
luminosity, Lx ~ L^oi. 



with luminosit}{3 as M ~ L^^^ ~ L^^^^, where the latter scaling uses a typical CAK 
power index a ^ 0.6. This then implies a super-linear scaling for X-ray to bolometric 
luminosity, Lj^ ~ Lj^'^p that is too steep to match the observed, near-linear law. 

In fact, the above scaling effectively assumes the shocks are radiative, with a cool- 
ing length that is much smaller than the local radius, ^ ^ r. In the opposite limit ^ » r, 
applicable to lower-density winds for which shocks cool by adiabatic expansion, the 
shock emission scales with the X-ray emission measure, EM ~ J p^dV, leading then 
to an even steeper scaling of X-ray vs. bolometric luminosity, ~ ~ L^-^^. 

Both these scalings ignore the effect of bound- free absorption of X-rays by the 
cool, u nshocked material that represents the bulk of the stellar wind. lOwocki & CohenI 
(119991 . OC99) showed that, in an EWS model in which the X-ray emission fraction 
drops with radius, accounting for wind absorption can lead to an observed X-ray lu- 
minosity that scales linearly with L^^*/- But while modern observations of spectrally 
resolved X-ray emission profiles by Chandra and XMM-Newton do indeed show the 
expected broadening from EWS, the relatively weak blue-red asymmetry indicates tha t 
such absorption effects are modest in even the densest winds (ICohen et al.ll2010l . l201lh . 
Since many stars following the L^-L^oi empirical law have weaker winds that are largely 
optically thin to X-rays, it now seems clear that absorption cannot explain this broad 
scaling. 



^For simplicity, we ignore here secondary scalings, e.g. of luminosity with mass, or wind speed with mass 
and radius. 
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The analysis here examines instead the role of radiative cooling, and associated 
thin-shell instabilities, in mixing shock-heated material with cooler gas, and thereby 
reducing and softening the overall X-ray emission. As summarized in figure 1, for a 
simple parameterization that this mixing reduction scales with a power (the "mixing 
exponent" m) of the cooling length, f^, we find that the linear Lj^-Lj^oi law can be 
reproduced by assuming m ^ 0.4. To lay the basis for deriving this result in §3, the next 
section (§2) first introduces a simple bridging form for emission between the radiative 
and adiabatic shock limits. 



2. Bridging Law for Adiabatic vs. Radiative Shock Emission 

Building upon the LDI-generated EWS scenario that is reviewed in this volume by 
Sundqvist et al., let us model the associated local X-ray emissivity within the wind as 

;7. = CpVv = Cp2.j-A_, (1) 

where p is the wind density, C is a constant that depends on the shock model and 
atomic physics, and /v represents a volume filling factor for X-ray emission. While 
previous work (e.g., OC99) has often directly parameterized this factor as following 
some specified radial function (e.g. power-law), the second equality in eqn. ([T]) casts 
this in terms of a "bridging law" between the limits for radiative <^ r) and adiabatic 
{i » r) shocks, where fq now represents some local "heating fraction", set by LDI- 
generated EWS. The cooling length itself scales as, 

rUc 1 

e= — = — . (2) 

P KcP 

where the cooling mass column depends on the energy of the EWS, and Kc = l/rric 
provides a convenient representation with units of opacity or mass-absorption coeffi- 
cient, e.g. cm^/g. In the simple model here, we assume that th e shock energy is fixed , 
and thus that Kc is spatially constant. From eqns. (18) and (22) of lAntokhin et al.l (|2004) . 
we find the numerical value 

Kc^ 190 E-lcm^/g, (3) 

where E^ev is the shock energy in keV. 

For X-rays emitted with photon energy comparable to the shock energy, the bound- 
free absorption opacity has roughly a similar inverse- square energy dependence, but 
with a numerical coefficient that is about a factor Kc/ki,/ ~ 8 smaller. For a wind 
with mass loss rate M and flow speed Voo, the transition from optically thick to thin 
X-ray emission can be characterized by the unit-optical-depth radius for bound-free 
absorption, 

R.^'-^. (4) 

In direct analogy, we can define a cliaracteristic adiabatic radius for transition from 
radiative to adiabatic cooling of the associated wind shocks, 

«,.^.140«e-^. (5) 
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where M_6 = M/10"^Mo/yr and Viooo = Voo/1000 km/s. For even the densest winds, 
the X-ray emission onset Ro ^ Ri, implying, as noted above, that wind absorption 
is at most a modest effect. But the stronger coefficient (Kc/ki,/ ~ 8) for radiative 
cooling means that such winds generally have Ra » Ro, implying that most O-star 
EWS remain radiative well above the wind acceleration region where they are gener- 
ated (IZhekov & Palla,2007.) . 



3. Thin-Shell Instabihty and Shock Mixing 

The inherent thinness of radiative sho ck cooling zones makes them subject to various 



thin-shell instabilities (IVishniadll994D . These can be expected to lead to an unknown, 
but potentially substantial, level of mixing between cool and hot material. Since cooler 
material radiates more efficiently, and in softer wavebands (toward the UV instead of 
X-rays), such mixing can significantly reduce the effective X-ray emission. While there 
have been some numerical simulations of the complex structure that arises from such 
instabilities (e.g. JWalder & Folinill 1998 1). there unfortunately does not yet appear to be 
any detailed study of how this can affect the net X-ray emission. 

To characterize the potential mixing effect on the Lj^-Lj^oi scaling, let us make the 
plausible ansatz that the reduction should, for shocks in the radiative limit ilr <^ 1, 
scale as some power of the cooling length ratio, £/r. To ensure that the mixing effect 
goes away in the adiabatic limit, we can (much as in eqn. [T]) assume a simple 'bridging 
law' scaling for a "mixing reduction factor" for X-rays, 

^ (1 + r/ir ' 

where the mixing exponent m > 0. To account for mixing within this model, we thus 
simply multiply the X-ray emissivity t]x in eqn. ([T]) by this mixing factor fjcm- 

As a specific, simple model, let us next also assume that, beyond some onset radius 
Ro, the X-ray heating fraction declines as power-law in radius, fq(r) = fqoiRolrY- Ne- 
glecting absorption, the X-ray luminosity can then be obtained from spherical volume 
integration of this X-ray emission. 



= 4;rC 



fxmfvP^r^dr = A7iCA- ; , (7) 

where Q = CfqoRl, and w(r) = V(r)/Voo is the scaled wind speed. For the standard 
13 = I velocity law, we have rw - r - R^, which even for general values of q and m 
allows analytic integration of (|7]) in terms of the Appell Hypergeometric function. As 
a specific example, for shock heating that declines with inverse radius {q = 1), direct 
numerical evaluation shows that the total X-ray luminosity is well approximated by a 
simple bridging law between the radiative and adiabatic limits. 
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M 

— I ; Ra<^Ro (9) 

— ; Ra^Ro, (10) 
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where the curly bracket term follows from straightforward integration of (7) for Ra 0, 
in which case the square bracket term just becomes unity; the latter scalings follow 
from limit evaluations of this square bracket, using the definition of Ra in eqn. (|5]). The 
transition Ra ~ Ro marks a kind of "sweet spot" for conversion of shock energy into X- 
ray emission; for lower density winds (Ra < Ro) much of that energy is lost to adiabatic 
expansion, while for higher density winds (Ra > Ro), it is lost to thin-shell mixing. 

In very dense winds with optically thick X-ray emission and so Ri > Ro, one 
can approximately account for the associated wind absorption through an "exospheric" 
approach (OC99) in which Ri simply replaces Ro as the lower bound for the integral in 
d?]), and thus also in (8). Since Ra/Ri = KclKbf ~ 8 » 1, the square bracket term just 
becomes a fixed constant, independent of M. Moreover, expansion of the curly bracket 
term now also makes the overall L^^ scaling independent of M for this fq ^ Ijr (q = \) 
emission case in the dense wind limit, Ra> R\^ Ro- 

4. Summary and Future Work 

The key result of this paper is that, in the common case of moderately dense winds 
with radiative shocks (Ra > Ro), thin-shell mixing can lead to a sub-linear scaling of 
the X-ray luminosity with the mass-loss rate, Ljc ~ (M/Voo)^'^- For a quite reasonable 
mixing exponent value m ^ I - a ^ 0.4, this then gives roughly the linear Lj^-Lj^oi law 
that is empirically observed for O-star X-rays. 

However we note that a similar mixing analysis could also be applied to model 
X-ray emission from colliding wind binaries, and their L^^ scaling with orbital separa- 
tion. Wide binaries with adiabatic shocks should still follow the usual inverse distance 
scaling, as directly confirmed by observations of multi-year-period elliptical systems 
like WR140 and rj Carianae (see review by Corcoran). But in close, short (day to week) 
period binaries with radiative shocks, mixing could reduce and limit the effective X-ray 
emission from the wind collision, and thus help explain why such systems often hardly 
exceed the L^^ ^ lO'^L^^/ scaling found for single stars (see review by Gagne). 

Finally, in addition to exploring such effects in colliding wind binaries, an over- 
riding priority for future work should be to carry out detailed simulations of the general 
effect of thin-shell mixing on X-ray emission, and specifically to examine the validity of 
this mixing-exponent ansatz for modeling the resulting scalings for X-ray luminosity. 
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